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Explicit expressions for the differential and total rates and emissivities of neutrino pairs from the 
photo-neutrino process e^ + 'y + u + i' in hot and dense matter are derived. Full information 

about the emitted neutrinos is retained by evaluating the squared matrix elements for this process 
which was hitherto bypassed through the use of Lenard's identity in obtaining the total neutrino 
emissivities. Accurate numerical results are presented for widely varying conditions of temperature 
and density. Analytical results helpful in understanding the qualitative behaviors of the rates and 
emissivities in limiting situations are derived. The corresponding production and absorption kernels 
in the source term of the Boltzmann equation for neutrino transport are developed. The appropriate 
Legendre coefficients of these kernels, in forms suitable for multigroup flux-limited diffusion schemes 
are also provided. 

PACS numbers: 52.27.Ep, 95.30.Cq, 12.15.Ji, 97.60. Bw 



I. INTRODUCTION 



In recent years, the study of neutrino emission, scattering, and absorption in matter at high density and/or tem- 
perature has gained prominence largely due to its inrportance in a wide range of astrophysical phenomena. Energy 
loss in degenerate helium cores of red giant stars P, Q , cooling in pre- white dwarf interiors Q , the short- and long- 
term cooling of neutron stars 0, , the deflagration stages of white dwarfs which may lead to type la supernovae 
0| , explosive stages of type II (core-collapse) supernovae , and thermal emission in accretion disks of gamma-ray 
bursters |^ [Tol | , are examples in which neutral and charged current weak interaction processes that involve neutrinos 
play a significant role. (The selected references contain more complete references to prior and ongoing work.) 

In unravelling the mechanism by which a type-II supernova explodes, the implementation of accurate neutrino 
transport has been realized to be critical ■ The basic microphysical inputs of accurate neutrino transport coupled 
in hydrodynamical situations are the differential neutrino production and absorption rates and their associated emis- 
sivities. The processes and precise forms in which such inputs are required for multienergy treatment of neutrinos for 
both sub-nuclear and super-nuclear densities (nuclear density po — 2.65 x 10^"^ g cm~'^) are detailed in Refs. 
At sub- nuclear densities, detailed differential information is available for pair production (e"*" + e~ —> v -\- v) |l2l| . 
nucleon bremsstrahlung (n-\-n^n-\-n-\-v-\-v) [T^ l. ^/-flavor production [vi + Di ^ i/j + Dj) p^ . and more recently 
for the plasma process (7* + D) [T^ . 

Our objectives in this work are to make available differential rates, emissivities, and the source and sink terms 
associated with the thermal photo production of neutrino-pairs in a Boltzmann transport treatment of the process 
6=*= + 7 ^ e"^ + V + v for which only the total rates and emissivities are available to date 0, 0, 0, |23, |^ l2^ . 
(It is important to note that in prior works, the energy and angular dependences of the emitted neutrinos were lost 
in simplifying the calculations of the total rates and emissivities; see Sec. II for details). In addition, we provide a 
qualitative physical understanding of the behavior of the neutrino emissivity for widely varying conditions of density 
and temperature. 

Section ^ is devoted to the derivation of working expressions for the differential rates and emissivities from the 
photo-neutrino process. The calculation of the hitherto unavailable squared matrix elements is outlined in Sec. II. A. 
Explicit expressions of these matrix elements, including those for the transverse and longitudinal components, are given 
in Appendix A. In Sec. II. B, the expressions for rates and emissivities are rendered in a form suitable for numerical 
calculations. The input photon dispersion relation is briefly discussed in Sec. II. C. Notes for obtaining accurate 
results from Monte Carlo integrations of the 8-dimcnsional integrals are provided in Sec. II. D and Appendix B where 
the resolution of numerical problems encountered in earlier works is addressed. Results of numerical calculations are 
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presented in Sec. III. The subsections here contain an analytical analysis of the qualitative behaviors of the total 
rates and emissivities for widely varying conditions of density and temperature. Sec. IIVI details the derivation of 
the production and absorption kernels in forms suitable for detailed calculations of neutrino transport along with 
numerical results for the leading Legendre coefficients. Sec. V contains a summary and a discussion of the relation of 
this work with those of prior works. Except when presenting numerical results, we use units in which h, c, and fcs 
are set to unity. 





FIG. 1: Leading order Feynman diagrams describing the emission of a neutrino pair from the photo process. The charged 
current W - exchange channel produces only Vei'e, whereas the neutral Z - exchange results in pairs of all three neutrino 
(e, n, and r) flavors. Contributions from positrons are obtained by the replacement e~ ^ e'^ . 



II. THE PHOTO-NEUTRINO PROCESS 

The leading order diagrams for the photoproduction of neutrino pairs, + 7 ^ + I'e (;j,t) + i^e (^i,t) 8i,rc shown 
in Fig. n The channel in which the exchange of the Z-boson occurs can produce any of the three species of 
neutrinos (I'c/j.r) and their anti-particles, whereas the channel in which the ly-boson is exchanged, only the Vei^e pair 
is produced. 

The total emissivity or the total energy carried away by the neutrino pair per unit volume per unit time from the 
photo-neutrino process is 



Q = 



(27r)3 2Ep J (27r)3 2lu J (27r)3 2Ep, J (27r)3 2Eg J {2tt)^ 2Eq, 
X {2i,f5^{j> + k-p' -q-q'){E, + E,,)-^Y.\^\''- (1) 



The first factor 2 accounts for the spin projections of the incoming electron. The factor ^ accounts for the polarizations 
of the incoming photon and the factor C arises from averaging over the spin projections of the outgoing electron and 
neutrinos. For the transverse polarization of the photon, ^ = 2 and = 4, while for the longitudinal polarization 
^ = 1 and C = 2. The index s keeps track of the spin projections of the initial and final electrons, while the index e 
does the same for the polarization states of the photon. The four-momenta of the participating particles are: 

p = {Ep,p) : incoming electron , 

p' = {Ep>,p'): outgoing electron , 

k = ((jj,k): incoming in — medium massive photon , 

q = (£'g,q) : outgoing neutrino , and 

q' = {Eq>,q'): outgoing anti — neutrino , (2) 
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where the energies and three-momenta are indicated in standard notation. In the entrance channel, electrons and 
photons are drawn from equilibrium Fermi-Dirac and Boson-Einstcin distribution functions 



expi ^1^1+1 



T 



, and 



cxp 



T 



c.-c.(k)-./|k| 



(3) 



respectively. The quantities /ie and T denote the electron chemical potential and temperature, while rrie and tUp stand 
for the electron mass and plasma frequency (this is related to the in-mcdium photon mass), respectively. The factor 
1 — Fe{Ep') accounts for the Pauli-blocking of the outgoing electrons. Blocking factors are absent for neutrinos, since 
the emitted (low-energy) neutrinos leave the production site without further interactions. 

The rate or number of neutrino pairs produced per unit volume per unit time, F, is given by an expression analogous 
to that in Eq. (Q, but without the factor Eg + Egi in the integrand. In prior works in which the total rates and 
emissivities were computed, the energy and angular dependences of the emitted neutrinos were eliminated by using 
Lenard's identity p^: 



(fq dPq' 

2Eq 2Eqi 



SHqt 



qW-^eiq^)i2q^q^ + q^gn, 



(4) 



where qt = q + q' = p + k — p' . Although the use of this identity simplifies considerably the calculation of the total 
emissivity, differential information about the neutrinos is entirely lost. On the other hand, calculations of differential 
rates and emissivities, such as 



dEq dEqi d{cos 6, 



and 



d^Q 



qq' ) 



dEq dEq' d{cOS 9qq> ) ' 



(5) 



where 9qqi is the angle between the neutrino pairs, entail the calculation of the relevant squared matrix element 
hitherto bypassed in obtaining the total rates and emissivities. We therefore turn now to the evaluation of the 
squared matrix element. 



A. Matrix Elements 



The Z— and W— exchange contributions to the total matrix element A4 = Mz +Mw can be combined by a Fierz 
transformation to yield 



leGp . I M 
M = — ^"e(p ) 



7"(Cv-Ca75) 



7 {(-^v - L-A75) 



2p - k + k- 



H + i- 



-2p' -k + k^ 



Ue{p) u^{q) 7a (1 - 75) Up {q') , 



where e is the charge of the electron and Gp is the weak (Fermi) coupling constant. For the + 7 
process, the numerical values of the vector and axial couplings, Cy and Ca, are 



(6) 



1 9 
Cv = - + 2 sin-^ Ow 



and Ca = - 
2 



where sin 9w = 0.226. The spin-summed squared matrix clement takes the form 



Tr if+me) 



+ me) 



Tr(^^7a(l-75)(/)(l +75)7/3) ,(7) 



where we have introduced the symbols 

Qi = p + fc , Q2 ^ p' - k , 

Pi = 2p-k + k^ and = -2p' ■ k + k^ . 

The emissivity from each of the -I- 7 — s- -f- i^{p,.T) + '^(p-.t) processes is obtained by the replacements 

Cv ^Cv-l and ^ Ca - 1 . 



(8) 



(9) 
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The total emissivity of all three neutrino flavors is obtained by the replacements 

^ + 2(Cy - 1)2 and ^ + 2{Ca - if ■ (10) 

For massless neutrinos, the trace over the outgoing neutrinos yields the familiar tensor 

Lap = 8 {qaq'p + q'^qf) ~ Qap q ■ q') + ie^a'^/S 9^9'" • (H) 

It remains then to contract this neutrino tensor with that obtained by performing the trace over terms that couple 
the electron with the photon. The result can be expressed as 



= 32e^Gl J2{^^v-CA)mlM^ + iCl + Cl) M+ + Cv CamA, 



(12) 



where the quantities A4-, M-\. and Mx depend on the scalar products of the various four momenta and the polar- 
ization of the photon. The remaining sum over the photon polarizations is performed in terms of its longitudinal and 
transverse components: 

V = -g^- + ^ - Pf" + Pr ■ (13) 

X=l 

The components of the transverse and longitudinal polarization tensors are 

p^. _ / for Aior!. = 

Pt - *,j- = 1,2,3 ' ^'^^ 

Pr = -a'-'+'^-Pr- (15) 

These polarization tensors satisfy the properties 

pMp p n pMP p pA* 



P^Plp. = -PL > PL = -1 , P^u - -2 • (16) 



[jl^''P r> pM pM 1 ^ . 

L ^Lpv — ^Ly ' ^Lp. ^ ^ ' ^Tp 

Explicitly, the transverse and the longitudinal components of the squared matrix element are 

Expressions for the quantities A^_, and A^x, and their transverse and longitudinal components are given in 
Appendix ^ 

B. Differential and Total Emissivities 

In this section, we derive expressions for the differential and total emissivities in forms that are suitable for numerical 
calculations. Similar quantities for the rates can be calculated analogously by dropping the factor Eq + Eqi in the 
integrand. We begin by rewriting the expression for the total emissivity in Eq. |^ as 



(2^)9 J 2E'qJ 2E^ J 2£p7 



Q = I ^ I ^ I ^ i^~Fe{E,,)] (Eq + Eq,) I{p\q,q') (18) 



I{p',q,q') = ji^l^^ I^F,{u;)F^{E,)SHp + k-P) ^l^''^'^^ (W) 

where the total four momentum and invariant squared mass are denoted by 

P = {E, P) = p + k ^ p' + q + q' and = . (20) 
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Note that the quantity I{p' , q, q') involves integrations over the incoming particles only. Utilizing the three-momentum 
delta function to integrate over the momentum p of the incoming electron, we obtain 



1 f°° IkP f^^ 

Iip\q,q')^-T^ ^d\k\ d^k d(cOS0fc)^7(^)^e(i?p)<5(i?-i?p-C.) 

H^T^) Jo ^p'^ Jo J-1 7^ 



(21) 



The energy delta function can be employed to perform integration over the angle 9k between P and k by using 



5{E-Ep-uo) _ 1 mi + \P\^ + \\^\' -{E-uo 



^ cos Ok - - ' ' ^ (22) 



which sets 



2Er, 2|P||k| V 2|P||k 



ml - AP -ml + 2Euj , , 

cos^fc = — ^TF^rn > (23) 



2|P||k 

where fc^ = m? = uj^ — |k|^. Integration over 0^ yields 



1 Ikl C'^'" 1 



(24) 



The condition | cos 6*^ | < 1 combined with Eq. (|23() can be used to establish the range in which either the energy or 
momentum of the incoming photon is able to conserve the total four-momentum P. However, the appropriate choice 

of the integration variable is dictated by the precise form of the dispersion relation uj = y^|k|^~i-m^. In the case that 
the mass of the photon is independent of its momentum, the momentum integration in Eq. H24() can be swapped with 
an energy integration using |k| d\\i\ ~ ujduj. The approximate dispersion relation, uj ~ \J |k|^ -I- uj^, for the transverse 
photon is an example of this situation in which mk = ujp. In this case, the photon energy must satisfy 

u' + Euj{ml^ ml - M^) + ^""^ ~ ~ + m^|P|^ < 0. (25) 

The roots of the above quadratic equation 



± _ E{AP + ml - ml) ± |P| + mi - mlY - AAPm,^ 



2 



2AP 



(26) 



specify the range in which the photon energy must lie. 

In the more general case that mk depends on the photon momentum |k|, it is convenient to retain momentum as 
the integration variable in Eq. H24|) . In this case, the range of momentum integration is obtained by the solutions of 



M^lkl^^^ - E{AP + ml - ml)\/ m"^^ + ml + E^ml + ^{Ap + ml - m^f = (27) 

which can be found by using an iterative procedure. 

In computing the differential and total emissivities or rates, it is desirable to allow the energies and the angle 
between the outgoing neutrinos to be as unrestricted as possible. The choice of coordinate axes that accomplishes 
this requirement for the outgoing particles and the concomitant restrictions on the incoming particles are described 
below. 

Outgoing Particles 

• The 3-momentum of one of the neutrinos is aligned along the z-axis. This fixes its 4-momentum to be 

q = (E,, q) = {E„ 0, 0, |q|) = (E,, 0, 0, E,) . (28) 

• The other neutrino is restricted to the x-z plane so that its 4-momentum is 

q' = {Eg', q') = {Eq>, |q'| sin 9qg, , 0, |q'| cos 9qq,) = {Eq' , Eq' sin 9qq, , 0, Eq' cos 9qq,) , (29) 

where 9qq, e [0, tt] is the angle between the two neutrinos. 
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• The energy Epi and the angles 9e e [0, tt], and 0e e [0, 27r] of the outgoing electron vary unrestricted in their 
domain so that 

p' = {Ep>, p') = {Ep', Ip'I sin de cos 0e , Ip'I sin 0^ sin , |p'| cos 6e) ■ (30) 

The specification of the outgoing momenta enables the determination of the total 4-momentum P, through Eq. I|20|) 
and hence |P| and angles 9p, (f>p can be obtained. 

Incoming Particles 

• Once P is known, the integration limits |k|^ for the photon momentum can be determined by using Eq. I|27|l . 
The simplest way to enforce these restrictions on |k|^ is to use a frame in which P is oriented along the z axis. 
For |k| to be in the allowed range in this frame, 9k is set from Eq. (|23|l and (j>k lies between and 27r. If energy 
is chosen to be the integration variable, Eq. (|26|l provides the appropriate range for the photon energy. 

• In order to use the newly determined incoming photon and electron momenta, we have to transform them back 
to the coordinate system determined by the outgoing neutrino-antineutrino pair. This is easily achieved by a 
rotation of k into the frame of the outgoing particles. The rotation is performed in two successive steps: first a 
positive rotation hy 9p about the x axis and then a subsequent rotation by 0p — 3tt/2 about the z-axis. 

• Finally, the 4-momentum of the incident electron p is found from p = P — k. 

The coordinate axes being chosen, the working expression for the total emissivity takes the form 

Q = J^J EqdE^j Eg, dE^, J cos 9g„ J l^d\p\' J ^dcos9e j dc^,[l - F,{Ep,)]{Eg + Eg,)I{p' ,q,q') , 

(31) 

where the quantity I{p' , q, q') is given in Eq. (|24|l . It is now straightforward to read off the differential emissivity from 
this result. Explicitly, 



d^Q 



dEgdEq, d{cos9gg,) (2 



2 roo ||-j'|2 /•! /■27r 

—E,E,, ^-^d\p\' dcos9j d<j>,[l~F,{Ep,)]{Eg+Eg,) I{p',q,q'). (32) 
'-'^) Jo Epi J_i Jq 



C. The Photon Dispersion Relation 



The dispersion relations of the photon in a plasma are commonly written as 

uj^{k) ^ k'^ +IlT{u}T{k),k) : Transverse 

u;l{k) = ^^I^ULiujLik),k) : Longitudinal, (33) 

where the transverse and longitudinal polarization functions 11^ and 11^ account for the effects of the medium on 
the photon. An extensive analysis of the exact dispersion relations for the densities and temperatures of interest 
here can be found in Refs. 0, Generally, calculations of 11^ and 11^ involve iterative procedures to solve 

either transcendental algebraic or integral equations. The transverse or longitudinal cmissivitics are then obtained by 



setting = y/IlT{k) or = y t^^f, — |kp, respectively, with the integration limits on the photon momentum |k|^ 
obtained iteratively from Eq. H27|) . This procedure, although straightforward, is time consuming. Fortunately, in the 
temperature and density ranges in which the photo-neutrino process dominates over the other competing processes, 
the leading order dispersion relations (see Fig. |2J) 



oj^ = uj'^ + |k|^ : Transverse (34) 

Luj^ = uip : Longitudinal (35) 

give an accurate representation of the full results. For the most part, therefore, we will present results by using these 
dispersion relations both because it is computationally faster and because it affords comparisons with results of earlier 
work. 
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FIG. 2: The leading order transverse and longitudinal photon dispersion relations ujt and lul used in the numerical evaluation 
of the neutrino emissivities. The momentum cut-ofFs {k^ /lJp)t,l depend on the density and temperature of the plasma. 



The relations in Eqs. and yield 



- |kr : 



Transverse 
Longitudinal . 



For the transverse photon, the momentum cutoffs are given by 



2AP 



1/2 



(36) 
(37) 



(38) 



where A = AP + uj^ — ml. 



For the longitudinal photon (plasmon), the momentum cutoffs are 



-,1/2 



(2|P|^ - 2L0pE + A)± w'(2|P|2 - 2ujpE + Af- - (4t^2^2 _ ^^^eA + A^) 



(39) 



The use of the leading order dispersion relations yields sufficiently accurate results except in regions where the 
plasma frequency becomes large so that the k integration enters the region where there are significant deviations from 
Eqs. jM | and p1|l . 

The transverse photon mass in Eq. H3()|) is momentum independent and therefore integration over the photon 
momentum can be replaced by integration over energy with limits provided by Eq. H26|l . We have verified that 
the two approaches yield the same result for this case. However, the longitudinal mass in Eq. H37() is momentum 
dependent and integrating over momentum is a better choice with the limits obtained iteratively from Eq. H27|) . 

The momentum cut-offs (fc^/wp)T,L are shown schematically in Fig. |21 Depending on the density and temperature 
of the plasma, the photo-neutrino rate and emissivity receive significant contributions from space-like momenta (i.e., 
|k| > Wp) in the longitudinal case. (In the nondegenerate regime, the time-like component is sub-dominant by more 
than two orders of magnitude.) This is in contrast to the decay of the longitudinal plasmon, 7* ^ v ~\- v, which is 
kinematically forbidden in space-like regions (i.e., |k| < ujp). 
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D. Notes For Numerical Integration 

For performing numerical integration, it is convenient to recast the integrals in Eq. ()31|1 into dimensionless forms 
by the variable transformations 

Eq=E^Xl, Eq, ^ X2, Ogg, ^ TT X3 , 

E,y = 4, + {E^, - 4,) X4, |k| - Ikf + (|k|+ - Ikf) xr , 

4>f, = 2tt 2:5, 9e = TT xg, 4)k = Stt (40) 

which results in 

/ 77^c\4 /•! /•! /•! /•! /•! 

Q = , "^o / 2:id.Ti / X2dx2 / sin6'qq, d.T3 / |p|' (E'^, - E'^,) da;4 / suvOedxQ 
(^■"■j Jo Jo Jo Jo Jo 

X / dx5 [1 - Ee(£;p')] (^9 + ^9') ^(J''' 9' (41) 
JO 

Iip',q,q') = ^ j^'M (|k|+-|kr)d.T7 j\xsF,iu;)F,{E,)± ^|A<^(^)p. (42) 

The 8-dimcnsional integral above is readily integrated by Monte Carlo methods. Although importance sampling 
reduces the variance of the final result, a flat sampling with suitable cut-offs of troublesome integrands yields equally 
good results. Such cut-offs can be easily identified through physical reasoning. 

In the case of degenerate electrons, the distribution function Fe{Ep, T) and the Pauli blocking factor I — Fe{Ep' , T) 
ensure that the relevant integrand peaks when the value of the electron energy is close to the chemical potential ^e- 
In other words, electrons primarily play the role of a spectator in this process, originating chiefly from the vicinity of 
the Fermi surface and reabsorbed into nearby available states. In the partially degenerate and nondegenerate cases, 
the width of the integrand is governed by the temperature T. We employ a flat sampling around the electron chemical 
potential /i^ and the width is taken to be ~ f OT, which gives the range of energy integration for the outgoing electron 
to be Ep = /ig ± f OT. Note, however, that rrie is the natural lower limit for the electron's energy, and should be 
employed when necessary. 

For neutrinos, the natural upper cutoff is given by the energy that is available in the medium, which in turn is 
given by the temperature. We can also expect that energy will be simply transformed from the photon to the two 
neutrinos, since the electron generally plays the role of a perturbed spectator. The energy released in the form of a 
I'D pair cannot greatly exceed the energy scale dictated by the temperature and setting E^ ~ (10 — I5)T makes a 
good choice for the integration cutoff. This can also confirmed by inspecting the integrand dQ/dE^: the integrand 
has a maximum around zero and decays roughly by 2 orders of magnitude for E^ > lOT. 

For certain physical conditions in the plasma, T > nie, and for densities in the range 10^" ^ ps^e/g cm~'^ ^ 10^^, 
numerical problems are encountered in using a Monte Carlo procedure to integrate Eq. H31|l (see also [TtL lisl HqI 
in which such difficulties have been reported). In Appendix B, we discuss the cause and remedy of this 
longstanding problem. 



III. RESULTS 



Our discussion will be restricted to the case of an equilibrium plasma in which the net negative electric charge 
of electrons and positrons is cancelled by a uniform positively charged background of protons, alpha particles, and 
heavier ions. The equation of state and the phase structure of matter, and the abundances of the various constituents 
including those of dripped neutrons at sub- nuclear densities are determined by the minimization of free energy [24| . 

We will present our results for the total and differential emissivities as a function of the mass density of protons in 
the plasma, psYe = mpUe, where mp is the proton mass. Ye = ng/uB is the net electron fraction {ns is the baryon 
number density) , and He is net electron number density 

1 f°° 

ne{T,fi,) = — dp {F,- - F,^) , (43) 
Jo 

which is simply the difference between the e~ and e+ number densities. (The quantity pbY,, is the same as p/ Pe used 
in prior works including Rcf. |25j.') Given n,., this expression can be inverted to determine the chemical potential /ie 
at a given temperature. 
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2 4 6 8 10 12 14 16 
logio(PBYe [g cm-3]) 

FIG. 3: The plasma frequency Up (upper panel) and the chemical potential /ie (lower panel) as functions of density for select 
temperatures. The vertical dashed lines in panel (a) indicate the densities at which o^p = T. The vertical dotted lines mark 
the turn-on densities at which the plasma frequency and the neutrino emissivity (Fig. 2J abruptly change their behavior from 
being independent of density. 

In order to gain a qualitative understanding of the basic features of Qt and Ql in terms of the intrinsic properties 
of the plasma, it is instructive to inspect the behaviors of the chemical potential //e and plasma frequency ujp (this is 
the characteristic energy scale generated by interactions in the medium) as T and psYe are varied. Figure El shows /ie 
(lower panel) obtained by inverting Eq. H43|l and ujp (upper panel) calculated from Eq. (59) of Ref. [13. The results 
of the total emissivities are readily interpreted on the basis of the trends observed in Fig. O 

Noteworthy features of the plasma frequency in panel (a) of Fig. |3|are: 

(1) LOp is independent of psYe tiU the turn-on density {pBYe)to is reached (this is at the root of why Qt and Ql are 
constant for psYe < (pBye)to), and 

(2) ujp shows a power-law increase for psYe > (pBi^e)to, the index depending both on the extent to which the 
plasma is in the non-degenerate, partially degenerate or degenerate regime and on whether electrons arc relativistic 
or nonrelativistic. 

The bold and light portions of the various curves in panel (b) of Fig. |21 mark the regions of densities for which 
/ie > T and p^ < T) , respectively. Inasmuch a,s pe ^ T indicates partial degeneracy of the plasma, the bold and light 
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transverse 
longitudinal 
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2 4 6 8 10 12 14 16 18 20 

FIG. 4: Individual contributions from the transverse and longitudinal channels to the neutrino emissivity as a function of 
baryon density at the indicated temperatures. The error bars show the variance of the Monte Carlo integration of Eq. 1411 . 
For densities in excess of nuclear density shown by the dotted vertical line, neutrino production from strongly interacting 
particles dominate over QED-plasma processes. 



portions refer to the degenerate and non-degenerate conditions, respectively. For reference, the electron mass, vifhich 
vifhen compared with or T determines the degree of relativity, is marked by the horizontal dotted line in this figure. 
The vertical dotted lines show the respective locations of the turn-on densities. 

The results of Monte Carlo integrations of the transverse and longitudinal total emissivities in Eq. I|41|l are shown 
as a function of baryon density in Fig. 0] for select temperatures. While the emissivities increase rapidly with 
temperature, their behavior with density is more intricate. At a fixed temperature, the important characteristics to 
note are: 

(1) Both Qt and Ql are independent of the density psYe until a turn-on density (pBYe)to is reached, 

(2) For densities larger than this turn-on density, Qt and Ql rise rapidly until they attain their peak values at a 
density (pBYe) pea^k, and 

(3) For psYe ^ (yOB^e)peak, the fall-off with density is exponential. 

A conspicuous feature to note is that with increasing temperature, the emissivities show a kinky behavior as they 
begin to approach their maximum values. 

In Fig. 13 the transverse and longitudinal rates, F^ and F^, are shown as functions of density and temperature. 
The qualitative features of the rates are similar to those noted above for the emissivities. 

The physical origins of the basic trends in the emissivities and rates are identified in the following section. 
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FIG. 5: Individual contributions from the transverse and longitudinal channels to the neutrino rates as a function of baryon 
density at the indicated temperatures. The error bars show the variance of the Monte Carlo integration. For densities in excess 
of nuclear density shown by the dotted vertical line, neutrino production from strongly interacting particles dominate over 
QED-plasma processes. 



A. Analytical Analysis of the Qualitative Behaviours 



In order to gain a qualitative, and in many cases quantitative, understanding of the neutrino emissivity as a function 
of temperature and density, it is useful to identify the dominant physical scales in the photoproduction of neutrino 
pairs under limiting situations. This can be achieved in the cases of (1) a degenerate plasma, which occurs for all 
temperatures at sufficiently high densities, and (2) a nondegenerate relativistic plasma in which the density is low, but 
the temperature is high. Finally, we turn to uncover the origin of the kinks that are conspicuous at high temperatures, 
but become less prominent at low temperatures. 

Since our goal here is to provide order of magnitude estimates to better understand the exact numerical results of 
the preceding section, we will simplify considerably the integrand in Eq. (|31|l by 

(1) estimating the dominant scales of energies and momenta of the incoming and outgoing particles, and 

(2) extensively use physically motivated approximations. 

These steps are helpful in taming the formidable 8-dimensional integral in Eq. H31|) which cannot be otherwise 
computed analytically. In the discussion that follows, we will focus on the dominant transverse channel; the analysis 
for the longitudinal channel can be carried along similar lines. 
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1. The Degenerate Case 
For all temperatures at sufficiently high net electron densities rie, the inequalities 

He > T, > rUe, 

ujp > T, tjp > nie (44) 

are satisfied. Thus, the electron chemical potential fi^ ~ (37r^ne) and the photon effective mass ojp ~ 4a/37r/ie ~ 
(e.g. see Refs. 0,113) are the dominant energy scales in the problem. Since T <^ /^e, the Pauli-blocking of 
outgoing electrons ensures that the participating electrons lie close to the Fermi surface. In other words, electrons 
are elastically scattered, exchanging only the 3-momentum with the photon and the outgoing neutrinos. Under such 
circumstances, we can expect the electron energies to be 

Ep^Ep, ^\p\^\p'\^^le (45) 

for both the incoming and outgoing electrons. This in turn implies that the photon transfers its entire energy to the 
vD pair. The photon energy and momentum are 

w ~ Wp ~ + Eq, and |k| ~ T . (46) 

The characteristic components of the outgoing neutrino 4-momentum then become Eg ~ |q| ~ LOpjl. These estimates 
for the relevant energies and momenta allow us to make an estimate of the squared matrix element. The generic form 
of the squared matrix element can be written as 

\M^f ~ 2567raG^(Cy + C^) x ^ — x {a 6*'^ order polynomial in 4 - momenta} , (47) 



(fc2 ± • fc)2 



where the dominant terms in the polynomial consist of 2 powers of fc and one power each of g, q', p, and -p' . The 
contribution from the term involving m1{Cy — C^) in the squared matrix element can be dropped, since rrie <SC jJ^e- 
By using power counting combined with our estimates of the relevant energies and momenta (q,; ~ Wp/2, pi ~ /^e, and 
k ~ Wp), the factor involving the polynomial can be approximated by ^ filujp/A. 
The factor 

]h " fc2 -2p'-k ^"^^^ 

needs more care, since for some ranges of temperature and density it can develop a resonant behavior (see Appendix 
B for a detailed analysis of the physical conditions for which this can occur). In the absence of such behavior, its 
magnitude can be estimated by averaging this factor over the appropriate angular variable. The largest contributions 
can be determined to be 

./f) ^ (m)^ ^^^^ 

Using these estimates, the squared transverse matrix element can be approximated by 



\M^\^ ~ WnaGliC^ + CI)ujI. (50) 



The remaining integrals in Eq. (|31|l for the total emissivity can be performed along the lines outlined in Sec. ^ The 
transverse emissivity then becomes 

4(.Z7r)» J Eqi Eq Ep, |F| J(_E_|P|)/2 

Since we expect that the most beneficial case for the emissivity is when Ep ~ Epi ~ /ig, we can set F{Ep) F{Epi). 
The dEp integration simply yields the range of Ep around /Xe for which the integrand is significant; i.e., / dEp ~ T. 
The factor 1/|P| = l/|p + k| can be approximated by 1//Ze, since |p| ~ fie ^ |k|. The dEp/ integrand is sharply 
peaked around Epi = He, so we can set Epi = jie everywhere except in the distribution functions, whence the dEpi 
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integration becomes straightforward. After introducing new variables x = Eq/T and y = Eqi/T, and the constant 
a = ujp/T, Eq. becomes 

Qt - // dxdyxy{x + y)e-^-+y^e{x + y~a), (52) 



(2^)« 

where the factor 9(x + y — a) arises from the condition Eq + Eqi = uj > uj.p. The integral above can be expressed as 

I{a) = J J dxdyxy{x + y)e-(^+y^e{x + y-a)^ -T{5,a), (53) 
where r(5, a) is the incomplete gamma function: 

Tin, a) = r dx = (a)""^ e"" ^ ^fl^ll^^ ' (54) 

The rightmost expression above is an asymptotic expansion. Since in this limit a = LUp/T ^ 1, we can keep only the 
leading order term in the asymptotic expansion to obtain 

n ~ 4 qCKC^ + Cj) 6 3 u;JT /rrx 

- 3 ^^0^ ^ ' ■ (^^^ 

Although the emissivity from this expression agrees very well with our numerical results at high densities, it should be 
emphasized that its utility lies in predicting qualitative trends, since the exact numerical prefactor depends strongly 
on the approximations employed. However, the Qt ^ ^^p^ q-'^p/t dependence can be employed to predict quanti- 
tatively the density at which the emissivity reaches its maximum value. From Eq. 1)55(1 . we obtain 



.i2( T ^ 



3 



(pB>;)max ^ 9.213 X 10^^ (^^j gcm-^ (56) 

At high temperatures (T > 10^ K) the emissivity has another, slightly higher maximum in the intermediate region 
where (/^e > T > ujp). This can significantly broaden the peak and mask the degenerate maximum. However, we 
can easily recognize the position of the peak corresponding to Eq. H5()|) by noting that a peak occurs just before the 
exponential decrease. These positions are accurately predicted by Eq. (|56|l with inputs for ujp and /ig from Fig. |21 At 
lower temperatures, or in the longitudinal case in the entire temperature range, this intermediate maximum is absent 
and we can clearly confirm the behaviour predicted by Eq. (|55|l . 

Employing the same physical reasoning and analytical approximations that were used in estimating the total 
transverse emissivity, the total transverse rate turns out to be 

^ 3 — — '^p^ n^^^p/T) 

3 (27r)6 ""P^ • ^^'> 



2. The Nondegenerate Case 

The nondegenerate situation occurs at sufficiently low densities for which fig — me ^ T. In this case, both Qt and 
Ql exhibit a plateau for temperatures T > 10^ K. This behavior of the emissivities is intimately connected with a 
similar behavior of the plasma frequency uip for the corresponding temperatures (see Fig. OJ. In this case, we can 
expect that a significant fraction of the electrons participate in energy exchange. 

For T > 10^" K, we can neglect the electron mass (me = 5.93 x 10^ K) in comparison to /ig. In the relativistic 
regime, net electron density and the plasma frequency are given by [T^ l25l | 



Me + J 



3^2 \^-<= . - •■ y - 3 

4a 
3^ 



2 4a/ 2 T^'^T'^\ 'iTTdT^ 
u^l = + (58) 
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where the rightmost expressions are vaHd for T 3> /ie- In this case, the characteristic electron energy and momentum 
are of order Ep ~ |p| ~ T, and the photon energy and momentum are w ~ |k| ~ T. 

Since the temperature becomes the dominant energy scale governing the neutrino-pair emission process, we can 
use dimensional analysis to establish that the neutrino emissivity and rate will scale with temperature as and T^, 
respectively. An analytical estimation of the numerical prefactors, however, needs some work. In this case, we have 
identified the leading term in the squared transverse matrix element to be 

> |7W I ~ 256TTaGp(Cy +C^) —— — -3 . (59) 

{2p'-k~k'^) 

The predominance of this term over all other terms is chiefly due to the resonant nature of the denominator. The 
angular dependences involving the outgoing neutrinos do not have a large effect, whence, to very good approximation 
we can write 

Y,\M^\' ^ 256 Tra G| [C'y + Cl) iuEpE.E,, ' ^\ . (60) 

^ {2p'-k-k^) 

We can further use the average value of 

{v'-k) 1 Ep,uj-\v'mx 



■ dx , (61) 



{2p' -k^k^f 2 7_i {2Ep,uj ^ 2\p'\\k\x - uj^y 

where x denotes the cosine of the angle between p' and k. The value of this integral does not depend sensitively on 
since it is negligibly small compared to the other terms in the denominator. We can therefore set it to zero without 
significantly changing the final result. In the relativistic regime, we can also ignore the electron mass compared to its 
typical momentum. We then find 

t^'^') . 1 1 .i„(i4).^, m 



{2p' ■ k - 2 AEp,uj \ujI J Ep-uj ' 

where in obtaining the rightmost result, we have utilized the average energy of an equilibrated nearly massless photon 
(w) = 2.7T and ujp ~ T/10 as appropriate for this regime. The working expression for the squared matrix element 
then becomes 

^lA^^I^ 2m.aGl{Cl + C\)^^^. (63) 



The identity 



'^'^ 6\q + g' - q,) E, E,, = ^ e{q^) {3qf - jqtp) , (64) 



2Eg 2Eq> " , V 24 



where the four- vector qt = {q^, qt) — q + q' = p + k — p' , helps us to write the emissivity and the rate as 
^Qt\ _ 32TT^aGj,{C^ + Cl) [ <Pp . d^p' d?k 



TtJ 3(27r)" 



f d-'p , , d^p' , , d-'k 



t2 , ,o> E„E„Er,' f E„ + u — E„ 



nB(u)) 

X Q{ql) {3qr - IqtH ^^^^ (^^' ' \ ) ■ (65) 

The remaining integrations are, however, complicated because of the complex integration boundaries arising from 

the Q{qt) factor. In order to simplify this restriction and the integration over the {Sq^^ — |cttP) factor, we use 
Q{qf) = 0(9t' — |qt|), and replace |qt| and |qtp by their average values: 

(|qt|> = I dy ^{El + El + 2E,E,,y) = + i|i q°j2 (66) 

(|qtP) = \j\y{El+El+2E,E,,v)^qf-2E,E,.^qf/2, (67) 
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where y denotes the cosine of the angle between the two outgoing neutrinos. In Eq. (|66|) . Ey and are the larger 
and smaller of the two energies Eg and Egi, respectively. In the last steps of Eqs. H66|) and H67|) . we have replaced the 
neutrino energies with their average value (?t /2. As a result, the following simplifications can be made: 

3qf - |qtP - {3qf - IqtP) = 5gf /2 . (68) 

We note that the condition 0(g°) = Q{Ep + w — Epi) has the physical interpretation that the outgoing electron energy 
cannot exceed the total incoming energy. 

The emissivity and rate can now be expressed in terms of the simple expressions 



Tt) 3(27r)6 {r^Xr 



where the dimcnsionlcss constants are 



/Xq\ r°° dxx^ 


f°° dyy 


f'+y dz 1 


Urj Jo e- + 1 j 


ey-ij 


e-^ + li 



{x + y~z) 



z)2 



(69) 



(70) 



and the integration variables are x = Ep/T, y = w/T, and z = Ep'/T. A numerical evaluation of these integrals 
yields 



Xq\ ^ /775.54 
Xt) ~ V136.50 

and the ratio Xq/Xt = 5.68. Analytical approximations to these integrals can be obtained by the replacement 

1 



(71) 



1 



1 (72) 



with the result 



= ^r(7)c(7)r(2)c(2) + |r(6)C(6)r(3)C(3) + |r(5)C(5)r(4)c(4) = 1002.1 (73) 

"^1 1 Q 7 

It = -r(6)C(6)r(2)C(2) + -r(5)C(5)r(3)C(3) + -r2(4)C2(4) = 172.74 (74) 

(r(n) and C,{n) are the Gamma and Riemann's zeta functions, respectively) and the ratio Xq/Xt — 5.8 which are 
reasonably close to the exact numerical results. Using the results in Eq. H71|l . the emissivity and the rate in the 
nondegenerate case are 

Qt ^ 4.35 X 10^2 (^^^ ergcm-3s-i, (75) 

/ rp \ 8 

-\27 I \ „_-3„-l 



Tt ^ 4.78 X 10^' (^^j cm-^s-\ (76) 

The above estimate for the emissivity is in good agreement with the exact numerical results presented in Fig. 0| 

As the density increases, the electron chemical potential and the plasma frequency both begin to increase, and 
become strongly density dependent. As a result, the emissivity also acquires a strong density dependence. We 
can expect this change of behaviour to occur when the term involving the chemical potential in Eq. H58|l becomes 
comparable to the term involving . In this case, 

(ne)to ~ 1^ or (77) 



( T \ 



3 



{psY,),^ 8 X 10^ j g cm-3 

, I ^^ _ { 10.7 at T = 8.6 MeV (10" K) 
logio(pBK,) - atT = 0.86MeV 



which agrees closely with the turn-on densities in Fig. 01 
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3. Intermediate Case 

In the case that fig > T > tOp, the emissivity Qt acquires a strong density dependence. In addition, a conspicuous 
maximum occurs prior to entering the strongly degenerate regime. In this density range, the dominant energies are 

electrons: E'p ~ |p| ~ /ig 

photon: w ~ |k| ~ T 
neutrinos: Eq~\(\\~T. (78) 

With these energy scales, the squared matrix element can be estimated to be 

Y1\M^\' - '256naG%iCl + Cl)J-^T^^^l (79) 

The phase space integration can be performed similarly to the degenerate case with the result 

" y W ^C(5)r(5,a;,/r). (80) 

In the region of maximum emissivity, the ratio LOp/T <C 1 for the two highest temperatures (10^" and 10^^ K). This 
allows us to expand the incomplete gamma function in terms of this small parameter so that 



768 aGliCl + Cj) TV" 
(2^)6 



For T = 8.6 MeV, the maximum occurs around psYe = 10^^ g cm^'^, where we have used /ie = 46.8 MeV and 
ujp = 2.74 MeV. By combining these values with Eq. (|81|l . the emissivity can be found to be logj^Q Qt = 32.06 in very 
good agreement that in Fig. 0] 

For T = 0.86 MeV, Qt exhibits a maximum around psYe = 10^ g cm^''. In this case, He = 4.71 MeV and 
u!p = 0.274 MeV. These values yield log^Q Qt = 23.0, which is slightly larger than the exact numerical result. 

The origin of the secondary peak lies in the resonant character of the factor 

^ = 1 ^.9^ 

for the case in which the plasma frequency becomes negligible, but Ep ~ |p| and w ~ |k|. Such conditions cannot 
be satisfied in a strongly degenerate medium in which w ~ 3> A: ~ T or at low temperatures for which Ep ~ 
and IpI ~ T ^ m^. Hence, as the temperature decreases this enhancement becomes reduced and at T < 10® K it is 
completely absent. 

This resonant structure also gives rise to numerical problems at high temperatures and is further discussed in 
Appendix Unless suitably accounted for, this factor enhances the variance of the Monte Carlo integration by large 
factors. 



B. Typical Neutrino Energies 

The mean neutrino plus anti-neutrino energy can be characterized by the ratio of the total emissivity to the total 
rate. In this section, we analyze 

{Eu,)t = ^ (83) 

for the transverse case. The analysis for the longitudinal case can be performed in a similar fashion. Utilizing the 
numerical results from Eq. (|41|l and its counterpart for the rate, the exact results for {E^p)t versus density at various 
temperatures are shown by the solid lines in Fig. The discussion below is aimed toward a qualitative understanding 
of the basic trends in this figure in limiting cases. 

In the degenerate case, the results in Eqs. (|55|l and (|57|) set the average neutrino pair energy to be 



{Evi,)rp ~ LOp 



(84) 
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in line with the assumption that Eq+E^i ~ LOp used in the approximation procedure. For the two highest temperatures 
r/Me V = 0.86 and 8.6 shown in Fig. El = 6.8r and IT, respectively. From the results in Eqs. (|7SJ) and lO 

m Sec. ITTT A 21 we obtain 



5.8r 



(85) 



in the non-degenerate relativistic case, which, considering the approximations made, is a reasonable estimate. 

It is instructive to compare the mean neutrino pair energy {Ei^p)j, with the characteristic energy of the in-medium 
photons in the plasma. The average photon energy in the plasma is given by 



[ljt) 



dk k'^ujT Fj{i^T,T) 



dk k^F^{ujT,T) 



8T 



(86) 



where the rightmost relation is obtained upon setting uj'^ 
for the extreme relativistic and nonrelativistic cases: 



Up (see also Ref. 0). Simpler results ensue 



(ujt) ^ 2.7 T for T > ujp 

3 

~ ujp + -T for T < Wp . 



(87) 
(88) 



The dashed curves in Fig. Elshow expectations based on Eq. (|86|l . 

Comparing the result in Eq. H85|l with that in Eq. (|84|) , we confirm that in the degenerate case the process proceeds 
as the decay of a massive photon into two neutrinos, with electrons and positrons emerging with modified angles in 
their final states. The numerical results in the figure also support this expectation. Notice, however, that as the 
density is progressively increased, the average neutrino energy becomes somewhat smaller than the average photon 
energy. This can be attributed to the fact that some of the available energy is taken by the outgoing electrons. 

In the relativistic case, the average neutrino pair energy (Eq. H85|l 'l is significantly enhanced relative to the average 
photon energy (Eq. H87|) '). This can be attributed to the different energy weightings of the squared matrix element in 
Qt and not accounted for in Eq. (|86|l and to the fact that electrons and positrons impart some of their energy to 
the neutrino pair. 



IV. KERNELS FOR NEUTRINO TRANSPORT CALCULATIONS 

The matrix elements derived in Sec. ^enable us to obtain the source and sink terms associated with the photo- 
neutrino process which are required in neutrino transport calculations. The evolution of the neutrino distribution 
function /, generally described by the Boltzmann transport equation in conjunction with hydrodynamical equations 
of motion together with baryon and lepton number conservation equations is 



dl 

dt 



df , difF^) 



BeaU) + BnesU) + BumU) + BtpU) . 



(89) 



Here, is the force acting on the particle and we have ignored general relativistic effects for simplicity (see, for 
example, Ref. ^3 for full details). The right hand side of the above equation is the neutrino source term in 
which, Bea^/) incorporates neutrino emission and absorption processes, B^Esif) accounts for the neutrino-electron 
scattering process, B^j^{f) includes scattering of neutrinos off nucleons and nuclei, and i?Tp(/) considers the thermal 
production and absorption of ncutrino-antineutrino pairs. 

In this section, we obtain the contribution from the photo- neutrino process to the thermal production term Btp- 
In doing so, we follow closely Ref. |l6l | where kernels for the plasma process and Ref. |l2l | in which neutrino pair 
production from e'^e^ annihilation were obtained. Suppressing the dependencies on (r, t) for notational simplicity, 
the source term for the photo-neutrino process can be written as 



B{f{flg,E,)) 



l-/(/.g,i?,) 

1 



1 f°° 



d^lq 



1 JO 



RP[Eq,Eq,, COS e, 



qql 



(2^ j / / ( Eq,Eq>, COS Oqq, ) f{Hq',Eq,) , 



l-J{iiq,,Eq,) 



(90) 
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FIG. 6: Average transverse photon and neutrino pair energies as functions of density and temperature. 

where the first and the second terms correspond to the source (neutrino gain) and sink (neutrino loss) terms, respec- 
tively. Angular variables ^9(9') = cos0q(q/) and ^^(q') are defined with respect to the z-axis that is locally set parallel 
to the outgoing radial vector r. The angle Bqqi between the neutrino and antineutrino pair is related to Qq and dqi 
through 

C01.Bqq, = flqflq, + \J 0- - " A*^') COs{(j)q ~ (j)q.) . (91) 

The production kernel is given by 

R^iEq,Eq,,coseqq,)=J J ^ J (2^)3 2E,, AEqE, ^ (^ + '-^ -'i-'i)-^Y:\M\ . 

(92) 

The corresponding expression for the absorption kernel R'^{Eq, Eq' , cos 6qq/) can be obtained by the replacements 

Fe{Ep) ^ 1 - Fe{Ep) , F^{uj) ^ 1 + F^{uj) , and 1 ~ Fe{Ep,) ^ Fe{Ep,) . 
The angular dependences in the kernels (Eq, Eqi , cos 9qq/) are often expressed in terms of Legendre polynomials as 

°° 2/ + 1 p 

R" (Eq, Eq,, cos 6qq,) = ^ __ $ » (i;„ ) ^'i (cos 0,,, ) , (93) 

(=0 

p 

where the Legendre coefficients (Eq, Eqi) depend exclusively on energies. 

From Eq. H92|l . it is evident that the kernels are related to the neutrino rates and emissivities. We first consider 
the production kernel RP {Eq, Eqi , cos 9qq,) . The corresponding analysis for the absorption kernel R"^{Eq, Eg' ,cos9qq,) 
can be made along the same lines. The neutrino production rate is given by 

2d^pFeiEp) fCd^kF^iiv) f d^p' [1 ~ Fe{Ep,)] f d^q 1 fd^q' 1 



(27r)3 2Ep J (27r)3 2uj J (27r)3 2Ep> J {2tt)^ 2Eq J (27r)3 2Eq^ 
2tt)^S\p^ 

d^q d^q' 



(2^)V(p + fc-p'-g-g')^El-^l'- (94) 



(27r)3 (27r)3 



RP iEq,Eq,, cos Oqq,), (95) 
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which defines the kernel RP{Eq, Eg, ,cos9qq,) and is to be identified with that in Eq. (|92|l . The emissivity Q can also 
be cast in terms of Rp using 

Q = I -^-^iEq+ Eq,)RP{Eq,Eq,,COSeqq,) . (96) 

Equations l|94(l and H96(l can be inverted to obtain 



RnEq,Eq,,COSeqq,) _ ^ dEqdE^d COS 9 qq, 

Stt^ d^Q 



(97) 



^2 ^2 (^^ ^ ) ^ cog ^^^^ ■ 

Combining Eq. (|32(l with the above equation, the production kernel can be written as 

RP {Eq,Eq,, COS 9qq,)^-—^-— ^dlpl'/ dcos9j d^^ - F,{Ep,)]I {p' , q , q' ) (98) 
[Zn)°EqEq, Jq Ep, J_i Jo 

where I{p',q,q') is given by 

1 /-oo ib-i /-s^ 1 

/b',9,g')-^(^y^ ^d|ki F,(.;) i^e(i?p) ^ Ei-^i'- (99) 

p 

The Legendre coefficients (i?g , iJg' ) are determined from 

<^[{Eq,Eq,) = J d{cOs9qq,)Pl{cOs9qq,)R"{Eq,Eq,,COs9qq,). (100) 

Using Eq. in Eq. (|1UI)|) . the Legendre coefficient for the production process can be expressed as 

3 pi foo I /|2 pi p2Tr 

^P{Eq,Eq,) ^ j—^—— d{cOS9qq,)Pl{cOS9qq,) ^rf|p|'/ d COS 9 , dcj, ,[\ - F,{Ep,)]I {p' , q, q' ) , 

{iTTfEqEq, J„i Jo ^p' J ~1 Jo 

(101) 

where I{p' , q, q') is given by Eq. H99|l . 

Numerical results 

The structure of the squared matrix element allows us to decompose the Legendre coefficients into parts that are 
symmetric and anti-symmetric in the energies Eq and Eq of the two outgoing neutrinos. For example, 

^^{Eq,Eq,) = {Cy + C\) X {part symmetric in [Eq, Eq,)} 
+{Cy — C\) X {part symmetric in {Eq, Eq,)} 

+CvCa X {part antisymmetric in {Eq, Eq,)} , (102) 

where only the parts symmetric in Eq and Eq contribute to the total emissivity Q and rate F, since the anti-symmetric 
part can be eliminated by the use of Lenard's identity fv\ IT^ I2H . Furthermore, it is clear from Eq. H93|) that only 
^Q{Eq, Eq,) contributes to the total emissivity Q and rate F; terms with / > 1 contribute only to the angular 
distribution. 

In order to ascertain the importance of the I > 1 coefficients relative to the / ~ coefficient, numerical integrations 
of the integrals in Eq. HlOlfl were performed for various values of I for the transverse case. Convergent numerical 
results were obtained by using three different methods: (1) a Monte Carlo procedure with uniform sampling using 
10^ — 10^ points, (2) the VEGAS Monte Carlo [23| procedure that employs stratified importance sampling, and (3) 
Gauss-Legendre quadrature with 8 and 16 points in each dimension. In regions of density and temperature where the 
resonant character of the squared matrix element does not become prominent (see Appendix B), the standard Monte 
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FIG. 7: Symmetric and anti-symmetric parts of the Legendre coefficients c]?^^''^'"' and cj)^''''"^™' in the production kernels for 
the transverse case. The neutrino energies Eg and Eqi are in MeV and the Legendre coefficients are in units of 10^^ erg~^ 
cm^ s~^. 



Carlo and VEGAS procedures yield similar variances (< 5%). However, a blind use of the VEGAS routine from 
standard libraries fails to evaluate the integrals properly in the high T and ps regions in which the resonant character 
of the squared matrix element strikes. As discussed in Appendix B, special care has to be exercised in treating this 
situation. 

Numerical results for the symmetric and anti-symmetric components of ^^{Eq, Eqi) and ^\{Eq, Eqi) are shown in 
Fig. |7|for the case T = 10^^ K= 8.62 MeV and psYe = 1 g cm~^. The results explicitly show the expected symmetry 
properties in Eq and Eqi . A comparison of the relative magnitudes in these four cases shows that <i>Q {^Eq , Eqi ) 
is the dominant term. The magnitude of ^^J^^^^^™) amounts to only 10% of the leading contribution. The 
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contributions of (ff* and ^J"-'^"^™-' are 6% and 3%, respectively. In physical terms, this means that neutrino-pair 
emission from the photo-neutrino process is dominantly isotropic. Therefore, depending on the required accuracy, 
(jjP(sym) jj^jgj-^^- |-,g adequate in practical applications. Note also that that the production kernels are negligible for 
energies Eq and Eq^ >10T. 



V. SUMMARY AND RELATION TO PREVIOUS WORKS 



In summary, we have calculated the differential rates and cmissivities of neutrino pairs from the photo-neutrino 
process e^-|-7 e^ + v+D in an equilibrium plasma for widely varying baryon densities and temperatures encountered 
in astrophysical systems. The new analytical expressions for the differential cmissivities yield total cmissivities that 
are consistent with those calculated in prior works [T^ Hol l20l I2lll2^ . In order to obtain these results, hitherto 
unavailable squared matrix elements for this process were calculated for the transverse, longitudinal, and mixed 
components. 

We have developed new analytical expressions for the total cmissivities in various limiting situations. These results 
help us to better understand qualitatively, and in many cases quantitatively, the scaling of the results with physical 
quantities such as the chemical potential and plasma frequency at each temperature and density. 

Using our results for the differential rates and cmissivities, we have calculated the production and absorption 
kernels in the source term of the Boltzmann equation employed in exact, albeit numerical, treatments of multienergy 
neutrino transport. We have also provided the appropriate Legendre coefficients of these kernels in forms suitable for 
multigroup flux-limited diffusion schemes. 

Beginning with the seminal work in Ref . [3| , all prior works 0, 0, |23, |^ have concentrated on providing 
improved total rates and cmissivities, but only for the transverse component. Our results for the total cmissivity 
show that the peak values of the longitudinal component arc nearly equal to those of the transverse component at all 
temperatures. (The mixed component vanishes exactly for the total cmissivity, but not for the differential cmissivity.) 
Through the years, prior works have persistently reported numerical difficulties in obtaining accurate results at high 
temperatures and high density. In this work, the cause of these difficulties has been traced to resonant factors that 
misbehave when Monte Carlo techniques are employed to perform multidimensional integrals. We have overcome these 
difficulties by utilizing a principal value prescription. The prinicpal new elements of this work are calculations of the 
full matrix elements, the neutrino differential rates and cmissivities, the analytical analysis of the total cmissivities 
in limiting situations, the development of production and absorption kernels in the source term of the Boltzmann 
equation for neutrino transport, and calculations of the appropriate Legendre coefficients of these kernels in forms 
suitable for multigroup flux-limited diffusion schemes. 
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APPENDIX A: SQUARED MATRIX ELEMENTS 



The squared matrix element in Eq. H12() contains the quantities A^_, and Ai^ which depend on the scalar 
products of the various four momenta and the polarization of the photon. The fruits of our toil arc summarized by 
the following explicit expressions for these quantities. 



M- = -q-q' 



PI 
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In writing the above expressions, we have employed the following notations: 

1 1 R _ 1 1 
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S{x,y,w,z) = (x-w) {yz) + {x-z) {y w) , 
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After summing over the polarization of the photon, the transverse and longitudinal components of the quantities 
^ T(L)^ j^^T{L) ^^^^ A^x"'"^^'' required in the calculation of the squared matrix elements in Eq. (|17|l are 
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The corresponding expressions for the fongitudinal component are 
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While summing over the photon polarizations, we have utilized the relations 



k2 

P^{x,y) = -x-y+ ^'^-'^l^-''^ -P^{x,y). (All) 
For completeness, the trace of six 7-matrices employed in the calculations above is given below: 



-^9ac,{9tii^9rl3 - gtirg^P + 9tifi9yr) 

+^9aT {9fiiy9<yl3 - 9f^a9^p + gt^pgya) ■ (A12) 



The trace with an additional 75 was computed by using 

) ■ (A13) 

APPENDIX B: TREATMENT OF RESONANT FACTORS 

In the h igh temperature regime, T > me, and for densities in the range 10^*^ PbYc/z cni^"^ ^ 10^^, prior 
works 0, 0, 0, 123, 0, 112 have reported numerical problems in calculating the neutrino emissivities through the 
use of Monte Carlo procedures. Our purpose here is to offer a solution to this longstanding problem. 

We have traced the origin of this problem to the factor 

1 ^ 1 ^ 1 

(32 k^~2p'-k P _2£;p'Cj + 2|k||p'|cos[<(k,p')] 

and its square in the squared matrix element. This factor gives rise to a sharp peak in the integrand for certain 
physical conditions in the plasma. In performing Monte Carlo integrations of Eq. H31|) , this feature causes a spurious 
increase in the emissivity as well as in its variance. 

At high temperatures, this numerical problem begins to occur at low densities, peaks at intermediate densities (in 
the case that ujp < T < ^g), and vanishes at sufficiently high densities. In what follows, we analyze the transverse 
component for which we use the lowest order dispersion relation. Generalization to the longitudinal component and/or 
the full dispersion relation is straightforward. The cases of T < mg and T > mg have distinctly different behaviors 
and are therefore considered separately. 

We begin by enquiring why this problem does not occur in the low temperature and low density regime (T <C rUe 
and fie ^ T). In this case, the energy of the outgoing electron is Ep/ ~ rrie and its momentum |p'| ~ 0. Furthermore, 
ujp ^ me (sec Fig. O. Therefore, we can drop the first and last terms in the denominator of Eq. (|B1|) so that 

^ ~ -2Ep'uj ^^^^ 

which is perfectly regular. 

With increasing density, the plasma begins to enter into the degenerate regime in which /ie ^ T, whereas fie ^ me. 
In this case, the energy and momentum of the outgoing electron can be approximated by Epi ~ |p'| ~ fie- For 
uJp 3> fie, significant contributions arise from to ujp and |k| ~ only, since photons with high energies/momenta are 
exponentially suppressed by the tail of the Bose- Einstein distribution function. In this regime, therefore, either |p'| 
or |k| is small. This allows us to drop the last term in the denominator of Eq. (|B1|) so that 

1 1 



(J2 i^p - 2fleUJp 

which is regular, since fie^ tOp- In physical terms, the fact that the energy of the outgoing electron is much larger 
than that of the photon ensures that Eq. (|Bip remains regular. 
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However, when the temperature exceeds the electron mass (T > me), Eq. (|B1|I can become extremely large. At 
low densities, which corresponds to the nondegenerate case (T ^ /ig), Up > rrif, (see Fig. O. In this regime, we can 
ignore neither the momentum of the electron compared to its energy nor the energy of the photon compared to that 
of the electron. The physical conditions are thus ripe for Eq. IjBip to blow up. The worst happens for intermediate 
densities for which uip <^ T <^ fj.^. In this case, we can ignore the plasma frequency, set Epi ~ |p'| ~ and w ~ |k| 
so that Eq. (|Bip takes the form 

^ ~ -2/^eC^(l + COs[<(k,p')]) ^^^^ 

which becomes large when the angle between k and p' approaches tt. 

As the density increases further, ^ |k| (see the all important Fig. |2|once more) and hence we can set lo = cUp for 
all relevant regions of the phase space not supressed by the Bose-Einstein distribution function. Since in this region, 
we have Epi ~ |p'| ~ ^ Wp, Eq. HB1|I turns regular again. 

In order to avoid the numerical problem associated with Eq. ljB4|) . it is necessary to tame the resonant character 
within the context of a Monte Carlo integration. Since we are evaluating a multidimensional integral, it is cumbersome 
to apply a different choice of variable sampling because the sampling function itself would depend on the remaining 
variables. The solution that is computationally simple is to evaluate the Cauchy principal value of the integral (see, 
for example, Ref. H^). This can be achieved simply by excluding the points when the angle between k and p' is 
within (a — e, a + e), where a is given by 

^Ep^oj - P 
2|k||p1 ■ 

In our computations, e ~ 0.01 proved to be effective in reducing the variance significantly. Note that this procedure 
can be adopted for any of the integration variables in Eq. |IB1|) , but at the expense of locating the lurking problems 
therein. 

We note, however that at high temperatures and high density, the plasmon and pair-neutrino neutrino processes 
generally dominate over the photo-neutrino process. 



A. V. Sweigart and P. G. Gross, Astrophys. Jl. Suppl. Ser. 36, 405 (1978). 
G. G. Raffelt, Phys. Rep. 333, 593 (2000). 

M. S. O'Brien and S. D. Kawaler, Astrophys. J. 539, 372 (2000). 

M. Prakash, J. M. Lattimer, R. F. Sawyer, and R. R. Volkas, Ann. Rev. Nucl. Part. Sci. 51, 295 (2001). 
D. G. Yakovlev, A. D. Kaminker, O. Y. Gnedin, and P. Haensel, Phys. Rep. 354, 1 (2001). 
W. Hillebrandt and J. C. Niemeyer, Ann. Rev. Astron. and Astrophys. 38, 191 (2000). 
K. Iwamoto, et al, Astrophys. Jl. Suppl. Ser. 125, 439 (1999). 
A. Burrows, Nature 403, 727 (2000). 

T. Di Matteo, R. Perna, and R. Narayan, Astrophys. J. 579, 706 (2002). 
K. Kohri and S. Mineshige, Astrophys. J. 577, 311 (2002). 

O. E. B. Messer, A. Mezzacappa, S. W. Bruenn, and M. W. Guidry, Astrophys. J. 507, 353 (1998); S. Yamada, H-Th 
Janka, H. Suzuki, Astron. Astrophys. 344, 533 (1999); A. Burrows, T. Young, P. Pinto, R. Eastman, T. A. Thompson, 
Astrophys. J. 539, 865 (2000); M. Rampp, H-Th Janka, Astrophys. J. 539, L33 (2000); M. Liebendoerfer, et al. Phys. 
Rev. D 63, 104003-1 (2001). 

S. W. Bruenn, Astrophys. Jl. Suppl. Se r. 58, 771 (19 85). 

A. Burrows and T. A. Thompson, ast ro-ph/0211404| 

S. Hannestadt and G. Raffeh, Astrophys. J. 507, 339 (1998). 

R. Buras, H-T. Janka, M. Th. Keil, G. G. Raffeh, and M. Rampp, Astrophys. J. 587, 320 (2003). 
S. Ratkovic, S. I. Dutta and M. Prakash, Phys. Rev. D 67, 123002 (2003). 
G. Beaudet, V. Petrosian and E. E. Salpeter, Phys. Rev. 154, 1445 (1966). 
G. Beaudet, V. Petrosian and E. E. Salpeter, Astrophys. J. 150, 979 (1967). 
D. A. Dicus, Phys. Rev. D 6, 941 (1972). 

J. R. Bond, Ph. D. Thesis: Neutrino production and transport during gravitational collapse (California Institute of Tech- 
nology) 1978. 

P. J. Schinder, D. N. Schramm, P. J. Wita, S. H. Margohs, D. L. Tubbs, Astrophys. J. 313, 531 (1987). 

N. Itoh, T. Adachi, M. Nakagawa, Y. Kohyama, and H. Munakata, Astrophys. J. 339, 354 (1989); and references therein 

for earlier works from the same group. 

A. Lenard, Physical Review, 90, 968, (1953). 

J. M. Lattimer and F. D. Swesty, Nucl. Phys. A 535, 331 (1991). 



26 



[25] E. Braaten, Phys. Rev. Lett. 66, 1655 (1991); E. Braaten and D. Segel, Phys. Rev. D 48, 1478 (1993). 
[26] P. J. Davis and P. Rabinowitz, Methods of Numerical Integration, 2"^ ed. (Academic Press, Inc., Orlando, USA, 1884). 
[27] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C: The Art of Scientific 
Computing, 2"'' edition, Cambridge University Press (1992). 



